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ABSTRACT 

We present a catalog of optical spectroscopic identifications of sources detected by Spitzer at 3.6 
or 24 /im down to -10 and ~280 ^Jy, respectively, in the SWIRE/ VMM- Aewton/ELAIS-Sl field and 
classified via line width analysis and diagnostic diagrams. A total of 1376 sources down to R ~ 24.2 
mag have been identified (1362 detected at 3.6 /xm, 419 at 24 /im, and 405 at both) by low-resolution 
optical spectroscopy carried out with FORS2, VIMOS, and EFOSC2 at the Very Large Telescope and 
3.6m ESO telescopes. The spectroscopic campaigns have been carried out over the central 0.6 deg 2 
area of ELAIS-S1 which, in particular, has also been observed by XMM-Newton and Chandra. We find 
the first direct optical spectroscopic evidence that the fraction of active galactic nuclei (AGN; mostly 
AGN2) increases with increasing F(2A fim)/F(R) ratio, reaching values of 70(±20)% in the range 
316<F(24 /xm)/F(i?)<1000. We present an Infrared Array Camera-Multiband Imaging Photometer 
color-color diagram able to separate AGN1 from obscured AGN2 candidates. After having corrected 
for the spectroscopic incompleteness of our sample, it results that the AGN fraction at F(24 /xm)~0.8 
mJy is ^22(±7)% and decreases slowly to ~19(±5)% down to F(24 /j,m)~0.3 mJy. 
Subject headings: galaxies: active - infrared: galaxies - surveys 



1. INTRODUCTION 

Multi- wavelength surveys are fundamental instruments 
to study the cosmological evolution of extragalactic 
sources. Each region of the electromagnetic spectrum 
provides information on different physical mechanisms 
which take place in the galaxies. The emission in the 
mid-infrared (MIR) band is mainly produced by dust 
heated either by the stars (particularly during periods of 
strong star formation activity) or by the emission orig- 
inating in an active galactic nucleus (AGN). In several 
cases, both causes contribute to the formation of the to- 
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tal spectral energy distribution (SED) in the MIR (see, 
e.g., Polletta et al. 2008; Donley et al. 2008 and refer- 
ences therein). 

In recent years, XMM-Newton- and C/iandra-based X- 
ray surveys have been able to delineate the density and 
evolution of X-ray unabsorbed and even moderately ab- 
sorbed AGN (Compton-thin; N H = 10 22 -10 24 cm" 2 ) up to 
redshift ^4 (see, e.g., Ueda et al. 2003; La Franca et al. 
2005). The fraction of absorbed AGN is observed to de- 
crease with increasing luminosity (e.g., Ueda et al. 2003; 
La Franca et al. 2005; Hasinger 2008; see also Lawrence 
&: Elvis 1982) and to increase with increasing redshift 
(La Franca et al. 2005; Ballantyne et al. 2006; Treister 
& Urry 2006; Hasinger 2008; Treister et al. 2009). In this 
framework, in order to reproduce the cosmic X-ray back- 
ground, a population of heavily X-ray absorbed AGN 
(Compton-thick; CT) as numerous as the Compton-thin 
one is expected. Because of the strong X-ray absorption, 
these CT AGNs are difficult to select even in the X-ray 
band. 

As the AGN (absorbed) activity is expected to con- 
tribute significantly to the MIR SED, observations in the 
MIR regime are potentially useful for selecting AGN, in- 
cluding those which are undetected in the X-ray band. 

The launch of the Spitzer telescope has allowed 
the construction of several large area multi-wavelength 
databases with MIR coverage, such as (for example) the 
9 deg 2 NOAO Deep Wide-Field Survey (NDWFS) of the 
Bootes field (Murray et al. 2005), the 3.7 deg 2 Extra- 
galactic First Look Survey (E-FLS; Fadda et al. 2006) 
and the smaller 0.5 deg 2 of the All- Wavelength Extended 
Groth Strip International Survey (AEGIS; Davis et al. 
2007). Indeed, the use of Spitzer colors has been re- 
vealed to be efficient (if not complete) in selecting AGN 
(e.g., Lacy et al. 2004, 2007; Stern et al. 2005). Us- 
ing both Spitzer and X-ray observations, Brand et al. 
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(2006), Treister et al. (2006), and Donley et al. (2008) 
have evaluated the fraction of AGN as a function of the 
24 /im flux, which results in being around 30%-45% at 
^3 mJy, decreasing to about 10% at ~0.4 mJy. 

Many studies have shown evidence of the existence of 
X-ray absorbed (Compton-thick) Spzteer-selected AGN 
(e.g., Donley et al. 2005, 2007; Alonso-Herrero et al. 
2006; Polletta et al. 2006; Steffen et al. 2007; Alexander 
et al. 2008). Maybe the most efficient criterion in select- 
ing X-ray-absorbed AGN is based on selecting sources 
with large F(24 /j,m)/F(R) ratios, in some cases linked 
to radio observations (e.g., Martmez-Sansigre et al. 2005, 
2007). Because of the faint optical fluxes these sources 
arc difficult to spectroscopically identify in the optical, 
and, in fact, the AGN classifications for this class of ob- 
jects are mainly based on MIR spectroscopy with the In- 
frared Spectrograph (IRS) on board Spitzer (e.g., Houck 
et al. 2005; Weedman et al. 2006; Yan et al 2007; Brand 
et al. 2008). Because of the low resolution of the IRS 
spectra, for many sources the classification and the red- 
shift measure are derived from the identification of the 
9.7 /im silicate absorption features on approximate power 
law spectral shapes. 

Using stacked Chandra images of sources with large 
F(24 nm)/F(R) ratios, Daddi et al. (2007) and Fiore 
et al. (2008, 2009) have shown that the average X-ray 
spectrum can be reproduced if a high percentage (~80%) 
of the sources in the sample are highly X-ray absorbed 
(even Compton-thick) AGN (but see Donley et al. (2008) 
and Pope et al. (2008) for partly different conclusions). 
This result is itself important for the understanding of 
AGN density evolution, but it needs to be complemented 
with a quantification of the real fraction of AGN among 
these sources. In fact, the above studies are mainly based 
on photometric redshifts and lack direct optical spectro- 
scopic identifications and classifications. Such measure- 
ments are crucial in order to understand how much of the 
average resulting hard X-ray spectrum is diluted and/or 
contaminated by starburst galaxies included in the sam- 
ples (see, e.g., Polletta et al. 2008; Donley et al. 2008). 

Brand et al. (2007) have studied in the near IR the 
spectrum of 10 ULIRG sources with 24 /im fluxes larger 
than 0.8 mJy and large ^(24 /j,m)/F(R) ratios and found 
that the SEDs were compatible with a mixed contribu- 
tion of AGN and starburst activity. Studying the SED 
and MIR spectrum of 21 obscured AGNs with large F(24 
[im)/F(R) ratios, and 24 /im fluxes larger than 1 mJy, 
Polletta et al. (2008) found that the contribution by the 
starburst component to the bolometric luminosity was 
below 20%. 

Dey et al. (2008) have measured spectroscopic red- 
shift for 86 very luminous dust-obscured galaxies with 
F(24 /im)>0.3 mJy and F(24 /im)/F(i?)>i000 in the 
Bootes field, and find a broad redshift distribution cen- 
tered at z ~ 2. Roughly half the redshifts are the results 
of Spitzer IRS observations, and half come from ground- 
based optical or NIR spectroscopy. However, no optical 
spectral classification of the galaxies has been discussed 
in their work. 

Here we present the optical spectroscopic identifica- 
tions and classification of 1376 sources of the central 0.6 
deg 2 of the ELAIS-S1 field included in the Spitzer Wide- 
area Infrared Extragalactic Survey (SWIRE; Lonsdale et 
al. 2003, 2004) and then use our database to estimate 



the fraction of AGN as a function of the MIR flux and 
the F(2A fim)/F(R) ratio. 

Throughout this work with "fraction of AGN", we 
mean the fraction of extragalactic sources in which is 
possible to reveal, via optical line width analysis and di- 
agnostic diagrams (see, e.g., Veilleux & Osterbrock 1987; 
Kewley et al. 2006), the presence of an active nucleus, 
regardless of its strength relative to the host galaxy. In 
the case of uncertain AGN signatures, we adopted a con- 
servative approach. Our database allows the first direct 
estimate (using optical spectra) of the fraction of AGN 
among the sources with large (>316) F(2A /im)/F(i?) 
ratios. 

We adopt a flat cosmology with Ho = 70 km s _1 
Mpc" 1 , ft M =0.30, and ft A =0.70. Magnitudes are given 
in the Vega system. Unless otherwise stated, uncertain- 
ties are quoted at the 68% (la) confidence level. 

2. MULTI- WAVELENGTH DATABASE 

The ELAIS-S1 field (center: a=00' l 35 m 00 ;s .0, S=- 
43°30'00"; J2000) is part of the SWIRE survey (Lons- 
dale et al. 2004) which is the largest Spitzer Legacy Pro- 
gram. SWIRE includes six high-latitude fields, totaling 
49 deg 2 (~7 of which on the ELAIS-S1 region) observed 
by Spitzer in four bands (3.6, 4.5, 5.8, and 8.0 /im) with 
the Infrared Array Camera (IRAC; Fazio et al. 2004), 
and in three bands (24, 70, and 160 /zm) with the Multi- 
band Imaging Photometer (MIPS; Rieke et al. 2004). 

Originally the ELAIS-S1 field was selected as one of 
the four fields of the European Large Area ISO Survey 
(ELAIS) covering a total of 12 deg 2 at 15 /im (Oliver 
et al. 2000; Rowan-Robinson et al. 2004). The size 
of the ELAIS-S1 field is 4 deg 2 and the 15 /im catalog 
published by Lari et al. (2001) contains 329 extragalactic 
sources over the flux range 0.5 — 100 mJy. The imaging 
and spectroscopic identification and classification of the 
15/im sources have been presented by La Franca et al. 
(2004; see also La Franca et al. 2007 and Gruppioni et 
al. 2008). 

An area of 3.9 deg 2 is covered by radio observa- 
tions obtained with Australian Telescope Compact Array 
(ATCA) down to Si^-ghz ^ 30 /zJy (la; Middelberg et 
al. 2008; see also Gruppioni et al. 1999). ELAIS-S1 
is also one of the targets selected by the GALEX Deep 
Survey, which has deeply covered the central part of SI 
in the far- and near-UV (Burgarella et al. 2005). 

2.1. Spitzer data 

The Spitzer/ SWIRE observations in ELAIS-S1 cover 
a ^7 deg 2 region larger than the whole 4 deg 2 ISO re- 
gion, in all the IRAC and MIPS bands, reaching ha sen- 
sitivities of 3.7, 5.4, 48, 37.8, 230, 18xl0 3 , and 150xl0 3 
/iJy in the 3.6, 4.5, 5.8, 8.0, 24, 70, and 160 /im chan- 
nels, respectively (Lonsdale et al. 2004). The SWIRE 
data in ELAIS-S1 have been released to the commu- 
nity in autumn 2005 (Data Release 3, DR3) through 
the query building GATOR at the NASA/Infrared Sci- 
ence Archive 16 . However, for the purposes of this work, 
we had access to the SWIRE working catalogs, which 
contain all of the sources in the public catalog, but 
reach deeper flux densities (on average about a factor 

16 http://irsa.ipac . caltech. edu/applicat ions/Gator/ 
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Fig. 1. — Position on sky of 3.6 /an (dots) and 24 /im (open 
squares) detected sources with a reliable redshift estimate. The 
circles indicate the area covered by the four XMM-Newton point- 
ings. The dashed and dotted lines show the statistically useful 
areas of the 24 fim and 3.6 /an spectroscopic samples, respectively 
(see Section 3.1). 

of 30% fainter). Details about the SWIRE data reduc- 
tion, that was carried out by the Spitzer Science Center 
and SWIRE team, can be found in the Data Release pa- 
per (Surace et al. 2005). Details on the generation of the 
Spitzer band merged catalog can be found in Gruppioni 
et al. (2008). 

While the original SWIRE/EL AIS-S1 field is -7 deg 2 , 
here we concentrate on a central rectangular area, of 
0.6 deg 2 , with limits 8.21°<a<9.245° and -43.91°<(5<- 
43.116° nearly including the four XMM-Newton observa- 
tions of the area described by Puccetti et al. (2006) (see 
Figure [T]). This area contains a total of 35021 sources 
detected by Spitzer, 32815 at 3.6 /mi and 2056 at 24 /im 
(1920 have been detected in both bands). Because of 
their bright flux limits, we excluded from the catalog the 
MIPS measures at 70 and 160 /im. 

2.2. Optical and NIR data 

B, V, and R images, down to B ~ 25, V ~ 25, 
and R ~ 24.5 (95% completeness) have been obtained 
with the WFI at the ESO 2.2m telescope within the 
framework of the ESO- Spitzer wide-area Imaging Sur- 
vey (ESIS; Berta et al. 2006; PI: Alberto Franceschini), 
while about 1 deg 2 of ELAIS-S1 has been covered by 
deep K and J bands exposures with SOFI at the ESO 
NTT telescope (I. Matute et al., in preparation). Re- 
cently, I- and z-band photometry, down to / ~ 23 and 
z ~ 22.5 (90% completeness), carried out with VIMOS 
at the Very Large Telescope (VLT), has been released 
(Berta et al. 2008). 

2.3. X-ray band data 

The first very shallow (~10~ 13 cgs; 2-10 keV) X-ray 
band observations on the ELAIS-S1 area were carried 
out with BeppoSAX and presented by Alexander et al. 



TABLE 1 
Number of good-quality 21 spectra 



Period 


Instrument 


Total 


3.6/im 


24/im 


V p_n a . 

A.&24/im 




15/iin Catal. b 


60 


54 


59 


8 


73-04 


VIMOS 


888 


885 


139 


35 


75-05 


VIMOS 


328 


327 


129 


18 


76-05 


EFOSC2 


20 


19 


12 


3 


77-06 


FORS2-SLIT 


11 


11 


11 


6 


78-06 


FORS2-MOS 


69 


66 


69 


13 


Total 




1376 


1362 


419 


83 



a We define as "good quality" those spectra with quality flag >1.5 
(sec scction l3.5t 

b La Franca et al. (2004, 2007) 



(2001). More recently, the central ~0.6 deg 2 region of 
ELAIS-S1 has been surveyed in the X-ray band with 
XMM-Newton (four pointings of about 70 "useful" ks 
each) with 478 sources detected, 395 in the soft (0.5 — 2 
keV) band down to a flux of 5.5 x 10~ 16 cgs and 205 in 
the hard (2 - 10 keV) band down to a flux of 2 x 10~ 15 
cgs (Puccetti et al. 2006). 

The regions with the highest XMM-Newton sensitivity 
(^65% of the full XMM-Newton area) were later target 
of 165 ks Chandra observations, reaching on-axis sensi- 
tivities of 2 x 10 -15 cgs (2-10 keV; S. Puccetti et al., in 
preparation), with the aim of obtaining precise positions 
for the X-ray sources. 

3. OPTICAL SPECTROSCOPY 

ELAIS-S1 was the target of several spectroscopic cam- 
paigns with ESO telescopes. The follow-up programs of 
the 15 /jm ELAIS-S1 sources provided 60 identifications 
of R<22 sources within our investigated area (La Franca 
etal. 2004,2007). In the period 2004-2006, five follow-up 
programs were accepted, with the aim of obtaining spec- 
troscopic identifications for the XMM-Newton, 3.6 and 
24 /im sources in the SWIRE/XMM-7Vewton/ELAIS-Sl 
area. The number of spectra obtained in each spectro- 
scopic run is shown in Table [TJ 

3.1. VIMOS observations 

Spectroscopic targets with a limiting magnitude of 
R ~ 24 were observed with the Visible MultiObject 
Spectrograph (VIMOS) at the VLT in multi-object spec- 
troscopy (MOS) mode, with the Low Resolution Red 
(LRR) grism (A/AA - 210), covering the 5500-9500 A 
wavelength range, with 1-4 hr exposure time per point- 
ing, totaling 66 hr of observing time. VIMOS data re- 
duction was carried out using the VIMOS Interactive 
Pipeline and Graphical Interface (VIPGI) developed by 
INAF Milano (Scodeggio et al. 2005). Two observing 
runs were carried out in the years 2004 and 2005. The 
first observing run was dedicated to identifying K-band 
and X-ray sources located in the western half of the area, 
as the SWIRE catalog was not yet available. The second 
run (as all the following runs) was dedicated to observing 
24 /im and X-ray sources, while the 3.6 /im sources were 
used as fillers of the masks. The observations were car- 
ried out in the eastern half of the area, but unfortunately 
it was not possible to observe the southeastern corner 
(see Figure[T]). In total 1473 spectra were collected: 1212 
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of good quality (quality flag >1.5; see section [3751 for a 
description of the quality flag). 

3.2. EFOSC2 and FORS2 (long slit) observations 

Four nights of visiting-mode observation in single-slit 
mode with EFOSC2 at the ESO-3.6m in La Silla allowed 
us to observe 20 optically bright (R<20) targets in the 
4000-9000 A wavelength range with grisms 6 and 13. 

Twenty hours of observing time was allocated with 
FORS2 at the VLT in long-slit mode. A total of 12 
good-quality spectra were obtained with the 150/ grism 
(A/AA ~ 260), covering the wavelength range 3700- 
10300 A. 

3.3. FORS2 (MOS) observations 

Spectroscopy of 69 faint sources has been carried out 
with FORS2 at the VLT in 2006 in MOS (MOvable-Slit) 
mode, with 2.5 hr exposure times per pointing, corre- 
sponding to a total of 28 hr of observing time. The same 
grism setup of the FORS2 Long-Slit observations was 
used (see the previous section). 

3.4. Data reduction 

The reduction process used standard MIDAS and 
IRAF facilities, except for VIMOS data which required 
a specific software package (VIPGI; Scodeggio et al. 
2005). The raw data were bias-subtracted, corrected for 
pixel-to-pixel variations (flat field), and eventually sky- 
subtracted. Wavelength calibrations were carried out by 
comparison with exposures of He and Ar lamps. Rela- 
tive flux calibration was carried out by observations of 
several spectrophotometric standard stars. 

3.5. Classification 

Spectroscopic redshifts were obtained using both the 
rvidlines IRAF tool and an instrument-optimized MI- 
DAS procedure, which allowed to compute the average 
of the redshifts corresponding to each measured line cen- 
ter. We assigned a quality flag to each redshift: [2] for 
the determinations based on at least two confirmed fea- 
tures (typically 3-4 or more absorption/emission lines 
besides the continuum shape), [1.5] for two very plau- 
sible features (besides the continuum shape), [1] for one 
clearly recognizable feature, [0.5] for tentative estimates 
of redshift. The spectroscopic catalogs we present in- 
clude only highly reliable redshift estimates, with quality 
flag > [1.5]. 

We classified the sources according to their optical 
spectral features in five broad classes. 

• Broad-line AGN (AGN1): type-1 AGN with broad 
(FWHM > 2000 km s -1 ) emission lines such as 
C iv A1549, C m] A1909, Mg n A2798, R/3, Ha. As 
a reference, in Figure O the composite spectrum 
of a sub-sample of 25 (l<z<3) 17 type-1 AGN is 
shown. 

• Narrow-line AGN (AGN2): type-2 AGN with 
narrow (FWHM < 2000 km s" 1 ) high-ionization 
emission lines (C iv A1549, C m] A1909, [Ne v] 

17 These redshift and wavelength intervals were chosen in order 
to maximize the number of overlapping spectra. 
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Fig. 2.— Composite spectra of 25 AGN1 (Kz<3), 7 high 
(1.6<z<2.6) and 12 middle (0.8<z<1.2) redshift AGN2, and 20 
ELG (0.6<z<1.2). 

A3426) or low-ionization ([O n] A3727, H/3, [O m] 
AA5007, 4959, Ha, [S a] A6725) emission lines 
with flux ratios indicating the presence of an AGN 
(e.g., Osterbrock 1989; Veilleux & Osterbrock 1987; 
Tresse et al. 1996). In Figure [2l the composite 
spectra of seven high (1.6<z<2.6) and 12 middle 
(0.8<z<1.2) redshift 17 AGN2 are shown. 

• Emission-line galaxies (ELG): sources with narrow 
emission lines, but no clear AGN signature in the 
optical spectra (see Figure [2]) . They show strong 
low-ionization emission lines that may be produced 
by thermal photons from hot stars. They often 
show also CaHK AA3933, 3969 absorption and the 
continuum break at 4000 A. This class is likely to 
include a fraction of NL AGN, possibly misclas- 
sified due to a small wavelength range (e.g., only 
[O n] A3727 visible), low signal-to-noise ratio (S/N) 
spectra, or dilution by the host galaxy emission. 
In Figure [H the composite spectrum of 20 middle 
(0.6<z<1.2) redshift 17 ELG is shown. 

• Absorption- line galaxies (GAL): sources with spec- 
tra typical of early- type galaxies, characterized by 
absorption features such as CaHK AA3933, 3969 
and strong 4000 A continuum break. 

• Stars: sources with stellar spectra (mainly of cold 
M-stars) , characterized by a blackbody continuum 
and absorption features. 
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TABLE 2 

The spectroscopic catalog 



IDesis 


R.A. 


Dccl. 


R 


z 


9 


els 




IDswiRE 


S 3.6fJ,m 


•524^1X1 


ID X 




[deg] 


[deg] 


[mag] 












[log(/LtJy)] 


[log(MJy)] 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


(9) 


(10) 


( n ) 


100^99 4Q 439M7 


8.84371 


43 491 ^9 


18.65 


0.225 


2 


4 


SWTRF3 


1003^99 ^0 439^17 3 


9 933^» 




311 


100^^91 ^zL 4'^9 c inQ M 

J UUOJ^ J. . 10^UU3 . O i- 




43 41 Q9k 


21.71 


0.636 


2 


2 


SWTRF3 


1003^91 3fi 439^00, 3 


1.7075 




301 




8 72240 


-43 30064 


20.65 


1.075 


2.0 


I 


SWTRF3 


1003453 39-431802 1 


1.8465 


2.5470 


258 


J003443. 36-431713. 39 


8.68062 


-43.28707 


18.42 


0.000 


2.0 


5 


SWIRE3. 


J003443.36-431713.3 


2.7933 




245 


J003407.39-430751.01 


8.53081 


-43.13084 


18.91 


0.188 


2.0 


2 


SWIRE3. 


J003407.41-430750.9 


2.2723 


2.1869 


185 


J003424.30-432037.21 


8.60127 


-43.34372 


21.47 


1.041 


2.0 


1 


SWIRE3. 


J003424.32-432037.2 


1.5117 


2.5663 


216 


J003429. 29-432409. 80 


8.62204 


-43.40269 


21.28 


1.065 


2.0 


1 


SWIRE3. 


J003429.30-432409.5 


1.7462 




229 


J003416. 03-433338. 06 


8.56669 


-43.56065 


22.17 


0.956 


2.0 


1 


SWIRE3. 


J003416.02-433338.2 


1.6669 


2.5285 


205 


J003320.65-433716.86 


8.33610 


-43.62132 


19.48 


0.287 


2.0 


2 


SWIRE3. 


J003320.68-433716.6 


2.1074 


2.4607 


80 


J003503. 93-432847.13 


8.76638 


-43.47975 


20.78 


1.108 


2.0 


1 


SWIRE3. 


.J003503.94-432846.9 


2.2380 


2.8999 


270 


J003716. 73-434151. 12 


9.31958 


-43.69757 


17.91 


0.226 


2.0 


3 


SWIRE3. 


.J003716.71-434151.1 


2.5961 


3.3699 
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Note. — Column 1: ESIS identification name from Bcrta ct al. (2006). Columns 2 and 3: coordinates (J2000). Column 4: i?-band Vega 
magnitude. Column 5: redshift. Column 6: redshift quality code: [2.0]— reliable, based on >2 confirmed lines; [1.5]— very plausible, based on 
two lines. Column 7: spectroscopic classification code: [1]— typc-1 AGN; [2]— type-2 AGN; [3]— ELG; [4]— normal galaxy; [5]— star. Column 8: 
SWIRE identification name of the most probable corresponding Spitzer source. Column 9: flux at 3.6 fim. Column 10: flux at 24 /im. Column 
11: XMM-Newton source name, sec Fcruglio et al. (2008). 

(This table is available in its entirety in a machine-readable form in the online edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content.) 



The wavelength range of our spectra (mostly in the 
interval 4500-9500 A), is wide enough to classify the 
sources. In the case of those narrow lines sources, where 
only the region of the [O n] A3727 line was visible, we 
used the presence (or not) of the high-ionization [Ne v] 
A3426 line to discriminate between AGN2 and ELG. In 
any case, the AGN2 classification has been conservatively 
assigned only in fully reliable conditions. Because of the 
presence of broad emission lines in their spectra, even 
in the case of low S/N data, we can exclude that type- 
1 AGN have been wrongly classified as either AGN2 or 
ELG. 

As far as the galaxies with no emission lines (passive) 
are concerned, we were not able to measure redshifts 
larger than one (with calcium-break redshifted out of the 
wavelength coverage). However, this bias should be small 
in the 24 /im sample, as the fraction of energy emitted 
in the MIR by the passive galaxies is smaller than in the 
other galaxies (see, e.g., the SED library shown in Fig- 
ure 1 of Polletta et al. 2007), and therefore the 24 /im 
sample should be scarcely populated by no emission line 
galaxies (see discussion in Section 4.1). This bias could 
instead affect the 3.6 /im sample 18 which, however, has 
not been scientifically used in this work. 

4. THE SPECTROSCOPIC CATALOG 

In Table [2j we present the whole spectroscopic cat- 
alog consisting of 1376 optical counterparts of 3.6 and 
24 fj,m Spitzer sources (whose redshift estimates have a 
quality flag > [1.5]). For every source we report: the 
ESIS catalog identification code (from Berta et al. 2006), 
sky coordinates in degree units (J2000) of the optical 
counterpart and its R-band magnitude, redshift, redshift 
quality code (only sources with 1.5 or 2.0 redshift quality 
are included), spectroscopic classification code ([l]=type- 
1 AGN; [2]=type-2 AGN; [3]=ELG; [4]=normal galaxy; 

18 Where indeed we found a fraction of classified no emission line 
galaxies four times larger than in the 24 (im sample; see discussion 
in Section 4.2. 
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Fig. 3. — 24 /an flux density distribution of all Spitzer sources 
(line on top), those with an i?-band detection (hatched region) and 
those with a reliable redshift estimate (darker shaded area). 

[5]=star). We report also the SWIRE identification code 
and the 3.6 and 24 flux densities, in log[^Jy] units, 
of the corresponding Spitzer sources. For those sources 
detected in the X-rays, the ID code from the catalog of 
Feruglio et al. (2008) is also reported. 

The spectroscopic sample of the (both Spitzer detected 
and not) XMM-Newton sources has been presented in a 
separate publication by Feruglio et al. (2008). There- 
fore, the spectroscopic identification of any X-ray de- 
tected Spitzer source appears in both catalogs. 

4.1. The 24 (J-rn sample 

Because of the incomplete execution of the second VI- 
MOS run, the statistically useful area of the field for 
spectroscopic identifications of the 24 /j,m sources is a 
rectangle with limits 8.21°<a<9.245° and -43.91°<(5<- 
43.116°, where all sources having a>9.0° and <5<-43.7° 
have been excluded (see Figure Q~|) . 

A total of 1932 24 /im sources were detected by Spitzer 
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Fig. 4. — Central: R magnitude as a function of the 24 fim flux. Those having a reliable redshift estimate are represented by symbols 
whose meaning is shown in the legend. Sources detected in the X-ray are represented by cyan crosses. The dashed lines indicate the 
loci with equal _F(24 fim)/F(R) ratio (whose values are reported in the labels). The vertical dotted line represents the 24 (im flux limit 
of 280 fim used to measure the fraction of AGN, while the horizontal dotted line represents the i?=24.2 mag limit of the spectroscopic 
identifications. Upper: fraction of the spectroscopically identified sources as a function of the 24 /an flux for the whole 24 /an sample 
(dashed line) and the i2-band detected 24 fim sample (continuous line). Right: fraction of the spectroscopically identified 24 fim sources 
as a function of the ij-band magnitude. 



within this 0.56 deg 2 studied area, 1512 of which (78%) 
have been identified in the R-band optical ESIS catalog. 
419 sources have been spectroscopically identified: 399 
result to be extragalactic sources and 20 stars. The ex- 
tragalactic sources were classified into 107 AGN (27.3%, 
52 AGN1 and 57 AGN2), 253 ELG (63.4%), and 37 GAL 
(9.3%). Figure [3] shows the histogram of the 24 /mi flux 
of all sources (top line), those with an R-band detection 
(hatched region) and those with a reliable redshift esti- 
mate (darker shaded area). 

In Figure |4j the distribution of the R magnitude as 
a function of the 24 /im flux is shown, with the indi- 
cation of the spectroscopically identified and the X-ray 
detected sources. The fraction of spectroscopically iden- 
tified sources decreases at both faint optical and 24 /tm 
fluxes. This trends should be taken into account when 
these data are used. In our analyses (see section[5|) , we se- 
lected a sub-sample with 24 //m fluxes brighter than 280 
/tJy and magnitudes brighter than R=24.2, which cor- 
responds to the largest possible sample where sufficient 
spectroscopic information exists (at these flux limits the 
fraction of spectroscopically identified sources is about 
10%). 

We have compared the spectroscopic sample (includ- 
ing only sources with a reliable redshift identification; see 
the previous section) with the "parent" 24 /im sample. 
In the upper left panel of Figure the fraction of the 24 
/im sources brighter than i?=24 mag which have been de- 



tected at 3.6 /im as a function of the 24 /im flux is shown. 
At 24 /tm fluxes brighter than 280 /tJy (once the uncer- 
tainties are taken into account), no significant difference 
is found between the spectroscopic sample and the par- 
ent 24 /tm sample: in both samples, in any bin more 
than 95% of the sources brighter than R=2A have been 
detected at 3.6 /tm. No significant difference is found, ei- 
ther, between the two samples, in the value of the F(24 
/im)/F(3.6 /xm) ratio as a function of the 24 /tm flux (see 
the lower left panel in Figure [5|) . We can then affirm 
that, although the fraction of spectroscopically identi- 
fied sources depend either from the optical and 24 /tm 
fluxes, the average MIR SED properties of the spectro- 
scopically identified sources are not significantly different 
from those of the parent 24 /tm sample (with R brighter 
than 24 mag). 

In the upper right panel of Figure El the i?-band distri- 
bution of the spectroscopically identified 24 /tm sources 
is shown. Unlike the identifications carried out with 
FORS2 and EFOSC2, where almost all spectra down to 
R ~ 24 mag are useful (quality-flag >1.5), at -R-band 
magnitudes fainter than R ~ 19 mag the fraction of low- 
quality VIMOS spectra increases with decreasing optical 
fluxes: it is about 20%, 35%, and 85% in the 21-22, 22- 
23, and 23-24 i?-band magnitude bins, respectively. As a 
consequence, most (~80%) of the spectroscopic identifi- 
cations at magnitudes fainter than i?=23 mag have been 
carried out only with FORS2, which is indeed much more 
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Fig. 5. — Upper left: fraction of the 24 (im sources brighter than 
i?=24 which have been detected at 3.6 fim as a function of the 24 
fim flux. The whole sample and the spectroscopically identified one 
are represented by open and filled circles, respectively. Lower left: 
ratio of the 24 fim and 3.6 fim flux densities as a function of the 
24 /im flux for the same samples. Symbols are as in the previous 
panel. Thick vertical error bars correspond to the (Iff) uncertain- 
ties on the mean, while the light vertical error bars show the (1<t) 
spread of the distributions. Upper right: R magnitude distribution 
of the spectroscopically observed 24/xm sources. The dotted line is 
the distribution of all the spectra, while the continuous line is the 
distribution of the "good-quality" (qual>1.5) spectra presented in 
this work. The shaded area is the histogram of those sources whose 
spectra have been mainly obtained with FORS2/VLT (some of the 
brightest with EFOSC2/3.6 m or DFOSC). Lower right: fraction 
of the "good-quality" (qual>1.5) spectra of the 24/im sources ac- 
cording to the spectroscopic class. Each histogram includes the 
classes below itself. 

efficient (because of better sky and fringe subtraction) 
at faint magnitudes. In the lower right panel of Figure 
the i?-band magnitude distribution of the fraction of 
each spectroscopic class among the good-quality spectra 
is shown. 

The redshift histogram of the spectroscopically iden- 
tified 24 fim sources and the 24 fim luminosity-redshift 
diagram are shown in Figure El The luminosities have 
been computed applying a fc-correction derived by a best 
SED fitting (F. La Franca et al., in preparation): follow- 
ing Fiore et al. (2008) we used a modified version of the 
SED library described in Polletta et al. (2007) and al- 
ready used for similar analyses by Lonsdale et al. (2004), 
Franceschini et al. (2005), Hatziminaoglou et al. (2005), 
Polletta et al. (2006), Weedman et al. (2006), and Tajer 
et al. (2007). The luminosity-redshift distribution is 
that typical of a flux-limited sample. On average, the 
most distant and luminous sources are the AGN1, then 
there are the AGN2, and finally the ELG. As this is also 
an optical limited sample (the spectroscopic identifica- 
tions reach R ~ 24.2 mag), the explanation of this be- 
havior requires a proper analysis (which is beyond the 
scope of this paper) of the bivariate MIR/optical lumi- 
nosity functions of each population, combined with the 
fc-corrections and all the selection effects. However, sev- 
eral works have been dedicated to the understanding of 
the shape of the 24 fim redshift distribution, discussing 
the combination of the effects of both the luminosity 
functions evolution and the MIR/optical SEDs (see, e.g., 
Perez- Gonzalez et al. 2005, Caputi et al. 2006, Desai 
et al. 2008). Desai et al. (2008) discuss in detail the 
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Fig. 6. — Top: redshift distribution of all 24 (im detected Spitzer 
sources. Those spectroscopically classified as AGN1, AGN2, and 
ELG are indicated in blue, red, and green colors, respectively. Bot- 
tom: \og(yh v ) 24 fim luminosity-redshift diagram of all spectro- 
scopically identified 24 fim sources. Blue circles, red squares, and 
green circles represent AGN1, AGN2, and ELG, respectively. 

redshift distribution of a sample of 591 24 /im sources, 
detected in the Bootes field of the NOAO Deep Wide- 
Field Survey, down to F(24 fim) = 300 fiJy and spec- 
troscopically identified down to R=25 mag (thus with 
flux limits similar to our sample) . They obtained optical 
spectroscopic redshift for 71% of the sources. Similarly 
to our sample, their redshift distribution shows a peak 
at z ~ 0.3 and a possible additional peak at z ~ 0.9 
(see their Figure 9). These peaks could be attributed 
to redshifted emission features in the SED which enter 
the 24 fim MIPS bandpass. According to the analysis of 
Desai et al. (2008), the z ~ 0.3 peak is difficult to be 
reproduced, as only the polycyclic aromatic hydrocarbon 
(PAH) emission features at 16.3 and 17 /im can be partly 
responsible of its presence, while the z ~ 0.9 peak could 
be attributed to the 12.7 fim PAH feature and the 12.8 
fim [Ne n] emission line passing trough the 24 fim band- 
pass. However, it is interesting to note that these two 
peaks are fairly well reproduced by the galaxy evolution 
model of Lagache et al. (2004) (see Figure 11 in Desai et 
al. 2008) where both normal and starburst galaxies are 
represented by luminosity-dependent SEDs. The model 
encompasses a strong rate of evolution of the luminosity 
density of starburst galaxies, peaking at z ~ 0.7 and re- 
maining constant up to z=A, while the normal galaxies 
evolve up to 2=0.4, after which their luminosity density 
remains constant. Analyzing a fraction of their sources 
without emission-line redshifts, Desai et al. (2008) find 
also weak evidence for another peak at z ~ 2 where about 
55% of the sources are AGN-dominated. 

At redshift larger than 1.0, our distribution is popu- 
lated by AGN. AGN1 show larger redshifts and lumi- 
nosities than AGN2: 80% of the AGN1 have redshift 
larger than ~0.9, while 80% of the AGN2 show redshift 
lower than ~1.1. This redshift distribution is not un- 
usual for an MIR sample combined to an optical flux 
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limit (introduced by the spectroscopic identification and 
classification process). See, e.g., the luminosity-rcdshift 
distribution of the 15 /im selected (F(15 /im)>500 /iJy) 
sample in the whole ELAIS-S1 region by Matute et al. 
(2006; their Figure 3) where the spectroscopic identifi- 
cations have been carried out down to R ~ 22 mag. In 
the 15 /jm sample used by Matute et al. (2006), 90% 
of the AGN1 show redshift larger than ~0.4, while 80% 
of the AGN2 have redshift lower than -0.4. The differ- 
ence between the redshift distribution of the two AGN 
populations is due to a combination of their luminosity 
functions with a difference in the average MIR/Optical 
ratio of their SED. Although both AGN1 and (with larger 
uncertainties) AGN2 show a strong luminosity evolution 
in the MIR (e.g., Matute et al. 2006 find a luminos- 
ity evolution at 15 /im of L(z) = L(0) x (1 + z) k with 
k ~ 2.9 for AGN1 and k ~ 1.8 - 2.6 for AGN2; see also 
Brown et al. (2006) for a recent estimate of the AGN1 
evolution at 24 /im), AGN2 have on average larger F(24 
fim)/F(R) ratios than AGN1 (see, e.g., La Franca et 
al. (2004) and the discussion in Section 5.1.1 of this 
paper). In our sample, about 40% of the AGN2 show 
F(2Afim)/F(R) ratios larger than 100, while only 15% 
of the AGN1 have F(24/im)/F(i?) ratios larger than this 
limit. A similar kind of difference between AGN1 and 
AGN2 is seen in the X-rays, as far as the X-ray to opti- 
cal ratio is concerned (see, e.g., Fiore et. al. 2003). This 
happens because, in AGN2, the optical AGN component 
is more often obscured, and the resulting optical spec- 
trum is the combination of the hosting galaxy emission 
with the AGN Narrow Line Region (when visible; see, 
e.g., Fiore et al. (2000, 2003), Cocchia et al. (2007), and 
Caccianiga et al. (2008) for a discussion on the X-ray 
Bright Optically Normal Galaxies-XBONG). The aver- 
age lower optical luminosity of AGN2 compared to the 
AGN1 (with the same 24 /im luminosity), combined to 
the optical R=24.2 mag spectroscopic limit, is therefore 
one of the most important agent for the average lower 
redshifts of the AGN2 sample (the luminosity function 
and the fc-corrections are the other most important ingre- 
dients). In summary, as the measure of the 24 /im AGN1 
and AGN2 luminosity function is beyond the scope of 
this paper, we can comment that the similarity between 
the AGN1 and AGN2 luminosity-redshift distribution of 
our sample and that observed at 15 /im (Matute et al. 
2006), suggests that it is qualitatively compatible with 
the previous estimates of the AGN evolution in the MIR 
(e.g., Matute et al. 2002, 2006; see also Brown et al. 
2006), where all these selection effects have been taken 
into account. 

4.2. The 3.6 fim sample 

Although the spectroscopic follow up was dedicated to 
the identification of the 24 /im and X-ray sources, a total 
of 1362 3.6 fim sources (881 not detected either in the 
24 /im or the X-ray bands) have been spectroscopically 
identified. These spectra come from the identifications of 
3.6 fim sources used as fillers of the VIMOS and FORS2 
masks, or from the identification of the isT-band (2.2 /im) 
sources of the western half area, observed during the first 
VIMOS run (see Section 3). We limit the following dis- 
cussion to the western half of the VIMOS area having 
limits 8.21°<a<8.74° and -43.91°<<5<-43.116° that can 
be treated as an unbiased representation of the parent 
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Fig. 7. — 3.6 (im flux distribution of all Spitzer sources (line on 
top), those with an i?-band detection (hatched region), and those 
with a reliable redshift estimate (darker shaded area). 

3.6 /mi sample and in which most of these sources (999) 
are indeed located (see Figure [1]). 

A total of 16782 3.6 /urn sources were detected by 
Spitzer within this 0.31 deg 2 area, 10794 of which (64%) 
have been identified in the -R-band optical ESIS catalog. 
999 sources have been spectroscopically identified: 762 
result to be extragalactic sources and 237 stars. The ex- 
tragalactic sources were classified into 111 AGN (14.6%, 
51 AGN1 and 60 AGN2), 358 ELG (47.0%), and 293 
GAL (38.5%). Figure [7] shows the 3.6 /im flux histogram 
of the sub-sample (line on top), those with an i?-band 
detection (hatched region) and those with a reliable red- 
shift estimate (darker shaded area). The fraction of 3.6 
/im sources with spectroscopic redshift is 22% for sources 
with F(3.6 /im)>50 /iJy and decreases at fainter fluxes 
(-1% at F(3.6 /tm)~30 /iJy). 

In Figure the distribution of the R magnitude as 
a function of the 3.6 /mi flux is shown, with the indi- 
cation of the spectroscopically identified and the X-ray 
detected sources. The fraction of spectroscopically iden- 
tified sources decreases at both faint optical and 3.6 /im 
fluxes. Unlike the 24 /im sample, a decrease in the frac- 
tion of identifications is evident at bright optical and 3.6 
/im fluxes too (see also Figure [7]). This is due to the 
fact that the spectroscopic identifications of the sources 
brighter than i?=18 were carried out by single slit obser- 
vations of the 24 /im sources only. 

This bias is evident in the upper left panel of Figure 
[9] where the fraction of the 3.6 /im sources brighter than 
R=24 which have been detected at 24 /im as a function 
of the 3.6 /jm flux is shown. In the 3.6 /im flux interval 
between 10 and 160 /iJy, there is no difference between 
the spectroscopic sample and the parent 3.6 /im sample. 
No difference is found either, in the same flux interval, 
between the ^(24 fxm) / F(3.6 /im) ratios of the two sam- 
ples (see the lower left panel in Figure [9]). On the other 
hand, at 3.6 /im fluxes brighter than 160 /iJy the fraction 
of 24 /im detected spectroscopically identified sources is 
significantly larger than the parent 3.6 /im sample. We 
can then conclude that, within the western half of the 
area, in each i?-band magnitude interval brighter than 
R — 24, the 3.6 /im spectroscopic catalog is a fair sample 
of the 3.6 /im sources only in the 10-160 /iJy flux range. 
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Fig. 8. — Central: R magnitude as a function of the 3.6 /im flux of the sources in the western half of the area. Those having a reliable 
redshift estimate are represented by symbols whose meaning is shown in the legend. Sources detected in the X-ray are indicated with 
cyan crosses. The two dotted vertical lines show the flux limits where, in each ij-band magnitude interval, the spectroscopic identifications 
are a fair ("gray") sample of the 3.6 ptm catalog (see Section 4.2). The horizontal dotted line represents the R=24.2 mag limit of the 
spectroscopic identifications. Upper: fraction of the spectroscopically identified sources as a function of the 3.6 /im flux for the whole 3.6 
Aim sample (dashed line) and the i2-band detected 3.6 Aim sample (continuous line). Right: fraction of the spectroscopically identified 3.6 
Atm sources as a function of the ij-band magnitude. 



In the upper right panel of Figure [9] the i?-band distri- 
bution of the spectroscopically identified 3.6 /im sources 
is shown, while in the lower right panel of Figure [9] the 
i?-band magnitude distribution of the fraction of each 
spectroscopic class among the good-quality spectra is 
shown. Unlike the spectroscopic 24 /jm sample (see Fig- 
ure [5] for comparison) even at faint i?-band magnitudes 
(R ~ 22 — 24), a large fraction (about 40%) of the sources 
are non-emission line galaxies (GAL) and stars. The red- 
shift histogram of the spectroscopically identified 3.6 /im 
sources is shown in Figure 1101 

4.2.1. IRAC colors 

With the launch of Spitzer, some color-color diagrams 
(mainly based on the IRAC photometry) have been pro- 
posed in order to quickly characterize the MIR SED 
properties of the sources and select the AGN (see, e.g., 
Lacy et al. 2004; Stern et al. 2005; Lacy et al. 2007). 
In Figure [TT1 the IRAC color-color diagram, as proposed 
by Lacy et al. (2004, 2007), of all the spectroscopically 
classified 3.6 sources is shown. The dashed line shows 
the original AGN selection criteria, calibrated on a Sloan 
Digital Sky Survey quasar sample and subsequently con- 
firmed with spectroscopic observations. The study of 
the relations among the average SED properties of the 
sources and their optical spectra is beyond the scope of 
this paper. However, we can note that, as in the above 
mentioned studies, also the sources of our sample occupy 



two branches, forming a "V" -shaped locus. Each spectro- 
scopic class preferentially occupies different regions of the 
diagram according to the energetically dominant compo- 
nent in the MIR regime (see, e.g., Sajina et al. 2005). In 
line with the above-mentioned studies, ~85% of the AGN 
(97% of the AGN1 and 70% of the AGN2 detected in all 
the four IRAC bands) in our sample are located inside 
the AGN selection region (the branch on the right) cor- 
responding to sources having red power-law MIR-SEDs 
whose slope increases moving to the upper right part of 
the diagram. ELG and GAL occupy (without solution of 
continuity) the upper and lower region, respectively, of 
the left branch. Stars are located at the bottom, almost 
in the center of the passive galaxy locus, as they have 
qualitatively similar SEDs (with an emission shortfall in 
the MIR). 

5. FRACTION OF AGN 

As discussed in the introduction, it is very interesting 
to measure the fraction of AGN among those sources 
showing large F(24 jiva)/F(R) ratios or as a function of 
the 24 /j,m flux. In order to use the most statistically 
useful sample, in the following analyses we will use the 
sub-sample of the spectroscopic catalog of the 24 /zm 
extragalactic sources corresponding to the intersection of 
the area of the spectroscopic follow up (see section 14. ip 
with the area covered by the XMM-Newton observations 
(see also Figure [T]). As discussed in section [4~Tl we select 
a sub-sample at 24 //m fluxes brighter than 280 fiJy and 
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Fig. 9. — Upper left: fraction of the 3.6 fim sources brighter than 
i?=24 which have been detected at 24 /im as a function of the 3.6 
jim flux. The whole sample and the spectroscopically identified one 
are represented by open and filled circles, respectively. Lower left: 
ratio of the 24 fim and 3.6 fim flux densities as a function of the 3.6 
fim flux for the same samples. Symbols are as in the previous panel. 
Thick vertical error bars correspond to the (lcr) uncertainties on 
the mean, while the light vertical error bars show the (lcr) spread 
of the distributions. Upper right: R magnitude distribution of 
the spectroscopically observed 3.6 fim sources. The dotted line is 
the distribution of all the spectra, while the continuous line is the 
distribution of the "good-quality" (qual>1.5) spectra presented in 
this work. The shaded area is the histogram of those sources whose 
spectra have been mainly obtained with FORS2/VLT (some of the 
brightest with EFOSC2/3.6 m or DFOSC). Lower right: fraction 
of the "good-quality" (qual>1.5) spectra of the 3.6 /mi sources 
according to the spectroscopic class. Each histogram includes the 
classes below itself. 
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Fig. 10. — Redshift distribution of the 3.6 fim detected Spitzer 
sources. Those spectroscopically classified as AGN1, AGN2, and 
ELG are indicated in blue, red, and green colors, respectively (as 
in Figure 6). 



magnitudes brighter than i?=24.2, which corresponds to 
the largest possible sample where sufficient spectroscopic 
information exists. These selections correspond to an 
area of 0.54 deg 2 . 




CW 300 

' . i .... i .... i .... i ... ' 
-0.5 0.5 1 

log(F 6B /F 36 ) 

Fig. 11. — IRAC color-color diagram of the spectroscopically 
classified 3.6 /an sources. Black points show all the sources de- 
tected in all the four IRAC bands. The dashed line shows the 
AGN selection criteria proposed by Lacy et al. (2004, 2007). Open 
black circles indicate those 2&fim sources having F(24 fim)/F(R) 
ratios larger than 300. For more information about the sources 
distribution, see, e.g., Sajina et al. (2005). 

5.1. Fraction of AGN as a function of the 
F(24 fim) /F(R) ratio 

The spectroscopic campaign has been carried out ob- 
serving sources detected by either XMM-Newton or 
Spitzer (or both). 83 of the 419 spectroscopically iden- 
tified 24 /iin sources are also included in the XMM- 
Newton detection catalog (see Table [lj . Therefore, in 
the estimate of the AGN fraction as a function of the 
F(24 fim)/F(R) ratio, a correction for the bias intro- 
duced by the inclusion of the X-ray sample has to be 
taken into account (indeed, as shown by Feruglio et al. 
(2008), most of the X-ray sources are AGN). If we call 
A^x (and iV XiBpec ) and A^NoX (and Af No x, sp cc) the number 
of Spitzer sources (and those spectroscopically identified) 
in a given MIR flux interval, detected and not detected 
in the X-ray band, respectively, the true fraction of AGN 
(Fr agn) among all (N^ot) the Spitzer sources is given by 

^ / iVx,AGN AT . A NoX ,AGN \ 

Ft AGN = -T7 N x + -r T A No x /N To t , 

\ J V X,spec JVNoX,spec / 

where iVx.AGN and A^noX.agn are the number of spec- 
troscopically identified AGN, detected and not detected 
in the X-ray band, respectively. 

We have verified that, although the fraction of spec- 
troscopic identifications decreases at faint i?-band mag- 
nitudes, in each F(24 fim)/ F(R) bin the fraction of spec- 
troscopic identifications is not significantly dependent on 
either the 24 fim flux or the i?-band magnitude (see Fig- 
urcll2|. Therefore, the spectroscopic identifications pro- 
vide a fair random sampling of the sources in each F(24 
fim)/F(R) bin. 

The fraction of AGN as a function of the F(24 
fim)/ F (R) ratio is shown in FigurefTSl while in Figure [HI 
the redshift of the spectroscopically identified sources as 
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Fig. 12. — Fraction of spectroscopically identified 24 /an sources 
as a function of the F(24 fim) flux divided into different bins of the 
F(24 fim)/ F(R) ratio (the dashed horizontal lines indicate mean 
values in each F(24 /j,m)/F(R) bin). Given the narrow dimension 
of the intervals in F(24 fim) / F(R), the overall constant trend in 
F(24 fim) is similar in F(R) as well. 
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Fig. 13.— Fraction of AGN as a function of the F(24 fim)/F(R) 
ratio (triangles). Filled circles represent the fraction of AGN1, 
while open squares represent the fraction of AGN2. Horizontal 
error bars represent the ratio interval; vertical bars represent the 
Poissonian errors propagated through the FYagn formula in Sec- 
tion 5.1. 



a function of the F(24 fim)/F(R) ratio is plotted. The 
fraction of AGN clearly increases with increasing F(24 
(jm)/F(R) ratio. While 85% of the AGN1 concentrates 
in the range l<log[i ;l (24 fim)/ F (R)]<2, which is typical 
of the local AGN1 SEDs (e.g., Spinoglio et al. 2002), 
the AGN2 fraction constantly increases with the F(24 
fim)/F(R) ratio: in the 2.5<log[ J F(24 fim)/F(R)]<3 in- 
terval (where a total of 18 sources have been spectroscop- 
ically identified) AGN2 outnumber AGN1 by a factor of 
~5, and the total fraction of AGN is 70(± 20)%. 

However, it should be noted that this statistics could 
be affected by some biases. In poor-quality data, emis- 
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Fig. 14. — Redshift of the spectroscopically identified sources as 
a function of the F(24 /im)/F(R) ratio. Filled circles and open 
squares represent AGN1 and AGN2, respectively. 

sion line spectra are more easily identified than no- 
emission line spectra. In the 2.5<\og[F (24 fim) /F(R)]<3 
interval, beside the 18 classified spectra, 8 further FORS2 
spectra (with similar i?-band distribution) missed an 
identification. These additional 8 spectra could be pop- 
ulated by the same mixture of AGN1, AGN2, ELG, and 
passive galaxies, but, under the most extreme assump- 
tion that all these objects were included among the non- 
AGN sample, the fraction of AGN in the 2.5<log[F (24 
fim)/F{R)}<3 interval would be 50%, instead of 70%. 
Therefore, although our estimate of a fraction of 70% 
AGN in the 2.5<log[F(24 fim)/F(R)]<3 interval has not 
to be considered an upper limit, it should be bear in 
mind that if this bias were present, the real dependence 
of the AGN fraction with increasing F(24 fim)/F(R) ra- 
tio would be weaker (but still present) than shown in 
Figure HH 

On the other hand, in Section 3.5, we have stated that a 
fraction (which is difficult to quantify) of the ELG classi- 
fied sources could instead reveal to be AGN2 if identified 
with higher S/N optical spectra. Although we are not 
able to quantify this bias, we can qualitatively expect 
it to be more relevant for the optically fainter sources, 
whose F(24 fim)/F(R) ratio is on average larger. In this 
case, the increase of the AGN fraction with increasing 
F(24 fim)/F(R) ratio would be even more relevant. 

In Figure [TBI the optical spectra of the 18 sources with 
2.5<log[F(24 fim)/F(R)]<3 are shown. The positions 
of typical AGN emission lines (even if not observed) are 
shown as a reference. Ten out of the 14 spectra with 
magnitudes larger than R—22.5 mag have been observed 
by FORS2. The AGN2 show larger [O m] A5007/H/3 and 
[O n] A3727/H/3 ratios than observed in the ELG, typical 
of AGN activity (e.g., Osterbrock 1989). Moreover, when 
observable, AGN2 show the high-ionization [Ne v] A3426 
emission line, which instead is not present in the ELG 
spectra (see also Figure [2]). 

The observed differences on the F(24 /j,m)/F(R) ratios 
of AGN1 and AGN2 is principally explained by the dif- 
ference between their optical luminosity. While in AGN1 
the optical luminosity originates directly in the nucleus, 
in AGN2, as the AGN is obscured, the optical luminosity 
comes from the lower luminous hosting galaxy. As shown 
in Figures [6] and [14] for non-AGNl sources the average 



12 



SACCHI ET AL. 



X 
> 

• r-H 
-1-1 

cd 
K 



4 


5 




4 


3 


5 




Q 
O 


2 


5 




2 


1 


5 




1 





5 


4 


5 




4 


3 


5 




Q 
O 


2 


5 




2 


1 


5 




1 





5 


4 


5 




4 


3 


5 




<3 
O 


2 


5 




2 


1 


5 




1 





5 


4 


5 




4 


3 


5 






2 


5 




2 


1 


5 




1 





5 


4 


5 




4 


3 


5 




3 


2 


5 




2 


1 


5 




1 





5 


4 


5 




4 


3 


5 




3 


2 


5 




2 


1 


5 




1 





5 




-, — , — , — , — I — , — , — , — I — I — I — , — ,- 

00039 R=21.2 z=0.861 AGN1 



00063 R=21.4 z = 1.176 AGN1 

is O Z 2 .£ 




-|^H 1 1 P — I 1 — f- 

14124 R=23.0 z=0.971 ELG 




frH 1 1 1 1 ' I ' I ' | ' — f — I 1 1- 



1$952 R=23.1 z=2.475 AGN2 T 18436 R=23.B z=1.025 AGN2 T 18842 R=23.0 z=1.069 ELG 

I 1 II § £ 



20568 R=23.6 z=0.972 ELG 



4-pi — i — I — | — i i ' i ' | ' — f — i — t^} 



22549 R=24.1 z=2.561 AGN2 



00044 R=21.6 z = 0.863 AGN2 



08036 R=24.0 z= 1.068 ELG 



I ' ' ' 



I I I I 



14943 R=23.3 z=1.119 AGN2 



f — I — I — I — '] — I — HH — | ' I ' I ' I — [ 



| 21571 R=22.6 z=1.076 AGN2 ! = 




-| — I — I — — I — V I ' | — I — I — I — |- 
22615 R=22.8 z = 1.238 AGN2 



j_L i i lL i i : i : . : . i l 



00062 R=21.8 z = 2.168 AGN1 



13348 R=22.9 z= 1.076 AGN2 



3 £, 



Z OSS £ £ 




I ' 1 1 1 I 



16373 R=22.5 z = 0.870 AGN2 




-| 1 1— I p-H p— I I^H f 



-| 1 1 h-^| V I ' | 1 1 1 1- 



22447 R=23.1 z= 1.353 ELG 
~ ^ r- < 



-| 1 1 1 1 ^ — I h-^ — i-Jh 1 1- 

22628 R=23.3 z=1.896 AGN 1 




I' i ■ . t I 



4000 6000 8000 



10" 40 6 8 1 4 4 6 8 1 0* 



Wavelength (A) 



Fig. 15. — Optical spectra of the 18 sources with 2.5<log(.F(24 /im)/ F(R))<3. The positions of typical AGN emission lines (even if not 
observed) are shown as a reference. 
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Fig. 16.— F(24 /j,m) /F(R) ratio as a function of the (IRAC) 
F(5.8 /im)/F(3.6 /im) ratio of 24 (im detected sources. The loci 
mostly populated by AGN1 and obscured AGN2 are delimited by 
dashed lines (the respective fractions are reported in the text). 
The larger F(24 fj,m)/F(R) ratio of AGN2 compared to AGN1 is 
due to the stronger obscuration in the optical for AGN2, while the 
F(5.8 /im)/F(3.6 jitm) separation between ELG and AGN is mainly 
due to the PAH features present in most IR galaxies combined 
with a redshift effect (for more information, see, e.g., Sajina et al. 
2005). The classification of the spectroscopically identified sources 
is shown in the legend. 

redshift increases with increasing F(24 /j,m)/F(R) ra- 
tio. All the sources with log [F (24 /im)/ F(R)]>2.5 have 
2>0.8, and vh v 24 /im luminosity larger than 10 44 erg 
s _1 . This behavior is due to the fact that those non- 
AGN1 sources showing larger F(24 fim)/F(R) ratio have 
preferentially the faintest optical magnitudes and thus 
the largest redshift. In contrast, AGN1 are more eas- 
ily found at large redshifts but with lower (roughly con- 
stant) F(24 nm)/F(R) ratios. In this case (as discussed 
in Section 4.1), this is caused by the evolution of AGN1 
luminosity functions (the density increases with increas- 
ing redshift 19 ) combined with the accessible volumes and 
^-corrections (see, e.g., Matute et al. 2002, 2006), while 
the average SED (and then the F(24 fim)/F(R) ratio) 
does not evolve significantly. 

5.1.1. IRAC-MIPS color-color selection of absorbed AGN 

As shown in Figure ITT1 the sources with large (>300) 
F(24 jj,m)/F{R) ratios are mostly located in the IRAC 
color-color region where AGNs are expected: i.e., in the 
branch where both the F(5.8 /im)/F(3.6 /im) and F(8.0 
/urn) / F (4.5 /nm) ratios are larger. However, it is evi- 
dent that the F(5.8 /im)/F(3.6 /im) ratio is more effi- 
cient than the F(8.0 /xm)/J 1 (4.5 /im) ratio in the AGN 
selection. It is therefore possible to create an optical- 
IRAC-MIPS color-color diagram, based on the F(5.8 
fj,m)/F(3.6 /^m) and F(24 fim)/F(R) ratios, which is 
highly efficient in selecting obscured AGN (AGN2). 

Figure [H shows how this optical-IRAC-MIPS dia- 

19 A luminosity evolution results in an increase of the density of 
sources with the same luminosity. 
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Fig. 17. — Fraction of optically classified AGN as a function of the 
24 ^m flux. The shaded area shows the lc confidence limits once 
the sources not detected in the R-band are included and assumed to 
contain a fraction of 68% of AGN, as measured in the 2.5<log[F(24 
fim)/ F(R)]<3 bin. Ovcrplotted are the estimates based on both X- 
ray and Spitzer data from Brand et al. (2006; squares) and Treister 
et al. (2006; triangles). The dashed line shows the expectations 
from the burst model of Pearson (2005). Horizontal error bars of 
our data represent the flux interval. 

gram is populated according to our spectroscopic clas- 
sification. The region with log[F(5.8 /im) / F(3.6 /xm)]>0 
and l<log[F(24 fim) / F (R)}<2 is mainly populated by 
AGN1: 65% are AGN1 and 12% are AGN2 (79% of 
all the AGN1 are located in this region). The re- 
gion with log[F(5.8 /rni)/F(3.6 fmx)]>0 and log[F(24 
lim)/F(R)}>2 is mainly populated by AGN2: 24% are 
AGN1 and 55% are AGN2 (36% of all the AGN2 are in- 
cluded in this area). Therefore, while the first region is 
quite efficient in selecting AGN1, the second one is use- 
ful in selecting obscured AGN2 with a moderate level of 
completeness. 

5.2. Fraction of AGN as a function of the 24 flux 

The measure of the fraction (and counts) of AGN as a 
function of the 24 /im flux is a very useful information to 
constrain the evolutionary models of AGN and starburst 
galaxies (see, e.g., Gruppioni et al. 2005; Franceschini et 
al. 2008). Unlike for the estimate of the AGN fraction as 
a function of the F(24 fj,m)/F(R) ratio (presented in the 
previous section), in this case, the presence in our sample 
of a decrease of the fraction of spectroscopically identified 
sources with decreasing i?-band magnitudes (see Figure 
[U and, on top of that, the lack of any identification at 
magnitudes fainter than i?=24.2, complicates the attain- 
ment of a reliable measure. 

In order to overcome these difficulties, and still obtain a 
sufficiently useful estimate, we have derived the fraction 
of AGN as a function of the 24 /im flux by assuming 
that it is mainly dependent on the ^(24 /j,m)/F(R) ratio 
as measured in the previous section. The AGN fraction 
was then derived by dividing each 24 /im flux interval 
into ^(24 /im)/F(i?) bins and multiplying the number 
of sources in each bin to its previously measured fraction 
of AGN (see Figure US]). 

In Figure 1171 the corresponding fraction of optically 
classified AGN as a function of the 24 /im flux is shown. 
The shaded area corresponds to the 1 a confidence lim- 
its once the sources not detected in the i?-band are also 
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included and assumed to contain a fraction of 68% of 
AGN, as measured for the objects with the largest F(24 
(jm)/F(R) ratio. The AGN fraction at F(24 /xm)~0.8 
mJy results to be ~22(±7)% and decreases slowly to 
~19(±5)% down to F(24 /W)~0.3 mJy. This estimate of 
the fraction of optically classified AGN, at 24 /im fluxes 
fainter than 0.8 mJy, is larger than the previous measures 
by Brand et al. (2006) and Treister et al. (2006) and the 
expectations of the burst model of Pearson (2005) (see 
Figure [T7)) . It should be noted that our measure of the 
AGN fraction is based on a direct spectroscopic opti- 
cal classification, while the estimates from Brand et al. 
(2006) are based on a classification which uses Spitzer 
colors (selecting sources in which the MIR emission is 
AGN dominated), and Treister et al. (2006) use a se- 
lection in the hard-X band and then corrects by the ab- 
sorbed AGN fraction expected to be missed assuming an 
Nh population distribution. The real fraction of AGN 
would be even larger if we took into account that a frac- 
tion (-10%) of X-ray (or even MIR) classified AGN do 
not show in the optical signatures of AGN activity (e.g., 
the XBONG; see Fiore et al. 2000, 2003; La Franca et 
al. 2007; Cocchia et al. 2007; Caccianiga et al. 2008; 
Gruppioni et al. 2008; Feruglio et al. 2008 and refer- 
ences therein). These objects usually show large values 
of both X-ray to optical and MIR to optical flux ratios. 

6. DISCUSSION AND CONCLUSIONS 

As discussed in the introduction, several studies have 
presented evidence that many galaxies showing SEDs 
with large F(2A /im)/F(R) ratios (e.g., Polletta et al. 
2008 use a threshold of F(24 fim) / F(R)> A00 and 24 //m 
fluxes larger than 1 mJy) harbor an AGN. These works 
were mainly based on low-resolution MIR spectroscopy 
or SED studies. X-ray stacking analyses of sources with 
F(24 /im)/F(R)>1000 have shown that their average 
X-ray spectrum is compatible with Compton-thin and 
Compton-thick absorbed sources (Daddi et al. 2007; 
Fiore et al. 2008, 2009), but the relevance of this re- 
sult depends on the fraction of AGN (and complemen- 
tary starburst galaxies) which contribute to the average 
X-ray spectrum (see Donley et al. 2008). 

Our catalog of 1376 optical spectroscopic identifica- 
tions of 3.6 and 24 /^m SWIRE sources, thanks to the 
use of FORS2 and VIMOS at VLT, has pushed the op- 
tical classification down to R — 24, and has allowed us 
to directly measure the fraction of AGN as a function of 
the F(24 fj,m)/F(R) ratio. We have then shown that, at 
24 /jm fluxes larger than 280 /iJy, in the range 316<-F(24 
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fim) /F{R) <1000, 70(±20)% of the sources show an op- 
tical AGN spectrum, and most of them are AGN2. In 
fact, the increase of the total fraction of AGN with the 
.F(24 /j,m)/F(R) ratio is caused by a strong increase of 
the fraction of AGN2, which populate more than 80% of 
the AGN found in the 316<F(24 /zm)/F(i?)<1000 range. 
This result is in agreement with the above-mentioned 
observed average X-ray absorbed spectrum reported for 
similar kind of sources. Indeed, in the framework of the 
classical unified AGN scenario X-ray absorbed spectra 
are typical of AGN2. 

At fainter optical luminosity (R>24) and then larger 
F(24 /j,m)/F(R) ratios the fraction of AGN might be 
even higher, but the optical spectroscopic identification 
of a statistical significant sample of these sources is dif- 
ficult to achieve with the currently available 8 m class 
telescopes. 

Our result confirms previous indications that the pop- 
ulation of sources showing large F(24 [im)/F(R) ratios 
are mainly AGN2 (e.g., Daddi et al. 2007; Fiore et al. 
2008, 2009), but as we were able to spectroscopically sep- 
arate the AGN(2) from the starburst galaxies, we are now 
able to measure how many and how much of these AGNs 
show X-ray absorption and, even more interestingly, if 
they include a significant fraction of the long searched 
and predicted Compton-thick AGN population. Such a 
kind of study requires the use of deep Chandra X-ray ob- 
servations which, in ELAIS-S1, partly already exist (S. 
Puccetti et al., in preparation), and will be the subject 
of a forthcoming paper. 
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